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^ ! Abstract 

. We report on a detailed analysis of generalization of the local adiabatic search 

I algorithm. Instead of evolving directly from an initial ground state Hamiltonian to 

' a solution Hamiltonian a different evolution path via a 'sure success' Hamiltonian is 

■ introduced and is shown that the time required to find an item in a database of size 
0^ ! N can be made to be almost independent of the size of the database, asymptotically 

I by modifying the Hamiltonian used to evolve the system. The complexity of the 

■ search algorithm in this and other similar approaches is shifted to construction of 
the Hamiltonian. 
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^ ' 1 Introduction 



One of the promising approaches to quantum computing is provided by the 
adiabatic evolution method introduced by Farhi et. al.[l] and elaborated and 
developed by several groups [1,2,3]. While the standard paradigm of quantum 
computation uses quantum gates (i.e., unitary operators) applied to quantum 
registers, the main idea of the adiabatic computing algorithms is to encode 
the initial state and solution state of the problem at hand in properly chosen 
Hamiltonians. Generally, the meaning of the "proper" is the requirement that 
the initial and the solution states are the ground states of an initial and 
final Hamiltonians. Then these two Hamiltonians are combined in a time- 
dependent way and system is let to evolve under the influence of this combined 
Hamiltonian. The speed of mixing is determined by the adiabatic theorem 
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of quantum mechanics [4,5]. According to this theorem, if the evolution of 
a quantum system is governed by a Hamiltonian that varies slowly enough, 
the system will stay near its instantaneous ground state. The general time- 
complexity of adiabatic quantum computation is still an open problem [6] . The 
complexity of the algorithm is determined by the time needed for switching 
from initial to solution state "slowly enough". The switching speed depends 
inversely on the size of the spectral gap between the instantaneous ground and 
first excited states of the Hamiltonian. 

One of the successful quantum algorithms is the database search algorithm of 
Grover [7], which has been examined in great details in many studies [8,9,10]. 
This algorithm shows a quadratic speedup in searching a database over clas- 
sical algorithms. Adiabatic search algorithm in its local adiabatic version was 
shown to provide similar speedup in searching, namely 0{-\/N) number of 
queries [3], where N is the number of items in the database. Also, some 
numerical evidence suggests that adiabatic algorithm might efficiently solve 
some instances of hard combinatorial problems, outperforming classical meth- 
ods [11,12,13,14]. 

Farhi et. al. [15] showed that one does not need to take a straight path from 
initial to final Hamiltonian. They showed that by applying a control Hamilton- 
ain at the evolution stage an algorithmic failure can be turned to a success for 
3-SAT problem. Recently, Boulatov and Smelyanskiy reported [16] a detailed 
theoretical study of using different paths for a version of 3-SAT problem. 

In the present article, we show that one can get a better than square-root 
speedup by using different path approach of Ref. [15]. We use part of 'sure 
success' Hamiltonian of Bae and Kwon [17] as the control Hamiltonian. Added 
control Hamiltonian keeps the minimum spectral gap a constant plus a system 

size dependent term which goes to zero as the system size tends to get close 
to infinity. It is shown that complexity of the proposed algorithm approaches 
a constant almost independent of the size of the system. 



2 Quantum Adiabatic Theorem 

Consider a quantum system in state \^{t)) which evolves according to the 
Schrodinger equation 

ij^\^{t))^H{t)\^{t)) 
where H{t) is the Hamiltonian governing the dynamics of the system (we 
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take ^ = 1). If the Hamiltonian is independent of time and the system is at its 
ground state, it will remain there forever. According to adiabatic theorem [4,5] 
if the Hamiltonian varies slowly enough, the state of the system will remain 
at the instantaneous ground state of the Hamiltonian for all t. 

Let \^k{t)) be the eigenstate of H{t), which satisfy 



H{t) \^k{t)) = Ek{t) \^k{t)) 

where Ek{t) is the corresponding eigenvalue and k labels the eigenvalues (A; = 
is the ground state) . 

The adiabacity conditions are related to the minimum energy gap of the spec- 
trum of H{t) and the transition to excited states of the system. These condi- 
tions can be made as follows: The minimum energy gap is defined as 



^min = ^min^ [Ei{t) - Eo{t)] 

and the maximum of the matrix element of dH/dt between the ground and 
the first excited state; 



D_ = max 



(1) 



Adiabatic theorem states that, if we prepare the system at time t = in its 
ground state |^o(0)) and let it evolve under the Hamiltonian H{t) for a time 
T, then 

\{^o{T)\^{T))f>l-e' 
provided that 



D 



max 



9rnin 



<e (2) 



where e -C 1. 



3 Local Adiabatic Evolution Algorithm 



The problem we consider here is to find a marked item in a set of items. 
Items in the system are labeled by n qubits, so the dimension of the Hilbert 
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space of the system is = 2". The initial state of the system is the uniform 
superposition of aU basis states, which also includes the unknown marked state 

\m) 



1 N-l 

iV'o) = ^ E (3) 

where \i) are the basis states with i — {),... ,N — 1. is the zero energy 
ground state of 



^o = /-|V'o)(V'o| 

whose ground state is |^o) with energy 0. Let the marked state Hamiltonian 
be 



I - \m){m\ 

whose ground state, the marked state |m), is unknown. The time-dependent 
Hamiltonian that interpolates between these two Hamiltonians is 



H{s) = {l-s)Ho + sHm (4) 

where the mixing parameter s = s{t) is a monotonic function with s(0) = and 
s(T) = 1, T being the required running time of the algorithm. The algorithm 
is to prepare the system in the state \ip{0)) = lipo) and then apply H{s) for a 
time interval T to end up with the state |m). Roland and Cerf [3] shown that 
complexity of this algorithm when applied in local adiabatic scheme scales 
with V^, similar to Grover's algorithm and this result is optimal for any 
chose of the time evolution of s{t). 

One can generalize the time- dependent Hamiltonian of Eq. 4 by adding to it 
a Hamiltonian of the form [18,17] 

Hd = a{s)\i/jo) {m\ + b{s)\m) {ipol 

where we have chosen a{s) — b{s) — s(l — s) so that the Hamiltonian at s = 
is Hq and at s = 1 it is Hm- Hd acts only during the adiabatic evolution of 
the system. So its effect is to change the path taken by the system for going 
from the initial state to final state \m). A similar extra Hamiltonian was 
shown to convert an algorithmic failure to a success by Farhi et.al. [15]. Now 
the time-dependent Hamiltonian of the system is 

H{s) = {l-s)Ho + sHrn-HD (5) 
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which in matrix representation can be written as 



' (s - \){x{s + - 1) (s - + x)\/\-x^^ 

^ s{l — s)x'^ 



[{s - l){s + x)Vl 



J 



where x = {m\i})Q) = 1/VN. Two lowest eigenvalues of this Hamiltonian are 
plotted as function of s in Fig. 1. These eigenvalues are separated by a gap 



g{s) = sjl - 4s(l - s) (1 - x{x + 1) - s(l - s)) (6) 
note that the minimum of g'(s) is at s = 1/2 with value g^min = 1/2 + 1/y/N. 
The eigenvectors corresponding to two lowest eigenvalues are 



l^o;s) 
\Ei;s)- 



{l-g{s)-2s)\m> +{s-l) (s + x) (l + vT^^) iV'o > 
^4 (s - if (s + xf (1 - x2) + (1 - g{s) - 2 s + 2 (s - 1) X (s + x)f 
(l + ^(s)-2s)|m> +(s-l) (s + x) (i + a/H^) IV'o > 
4 (s - 1)^ (s + a;)^ (1 - x"^) + (1 + ^(s) - 2 s + 2 (s - 1) a; (s + x))^ 



1 .0 ^ 
0.8 : 

0.6 : 

0.4 - 



0.2 




Fig. 1. Eigenvalues of the two lowest states of the time-dependent Hamiltonian H{s) 
function of reduced time s for N = 64. 

The matrix element in Eq. 1 can be written as 



'dH\ _ ds ldH\ 
. dt 1 1,0 dt\ds / ^ Q 



ds 



(7) 
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(//„-//o + (1-2s)//d>i,o<1 [3]. 

In the the local adiabatic approximation one tries to adjust the evolution rate 
of the s for infinitesimal time intervals dt and apply the adiabacity condi- 
tion locally to each of these intervals. The adiabacity condition for the local 
adiabatic evolution would be 



ds 
'dt 



< € 



As) 



/dH\ 
\ /l,0 



for all times t. Using Eq. 6, the evolution of the Hamiltonian is chosen to 
evolve at a rate that is solution of 



g.es^W^e 1-4.(1-.) + 



(8) 



Integral of this equation involves the logs of the roots of g'^{s) and is displayed 
in Eq. 9. s{t) obtained from this integral is shown in Fig. 2 along with the s{t) 
of Roland and Cerf 's algorithm. 



V2h{N) \k2{N) 



tanh ^ 



2V2k2iN)s 
2-k2{N)'^-As 



ki{N) 



tanh ^ 



2v^A:i(iV)s 
2-ki{N)^-As 



where 



h{N) = -^/(l-N)(2VN+l) 



M^) = ^-2-4MiV) + A + -^ 
k2{N) = l-2 + ih{N) + ^+ ^ 



N y/N 



As can be seen from the figure, added term Hd changes H{s) much faster 
compared to the original Hamiltonian of Ref. [3]. The computation time of 
the algorithm is evaluated by taking s = 1 in the solution of Eq. 8, and is 
plotted against in Fig. 3 along with the same quantity from Ref. [3]. The 
asymptotic form of Eq. 8 for s = 1 is given by 



T 



1 + 



TT 



+ o (i/^/n) 



(10) 
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Fig. 2. Dynamic evolution of the Hamiltonian that drives the initial ground state to 
the solution state: plot of the evolution function s{t) for N = 64. The dotted line is 
for the same quantity calculated by using the Hamiltonian of Ref . [3] . 



which show that as the system size approaches infinity, the required running 
time approaches the constant (l + j^. 
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Fig. 3. Computation time vs the size of the problem. The dashed and dotted line is 
for the \/iV and Hamiltonian of Ref. [3] , respectively. 



The main factor that determines the required evolution time T is the behavior 
of minimum gap g-aan- This energy difference usually reaches its minimum at 
an avoided level-crossing. If this approach is exponentially fast with a growing 
particle number then the algorithm is not efficient. The effect of extra term 
Hd is to keep these two levels far apart. Our result seem to be similar to one 
obtained by Das et. al. [2] who use Eq. 4 with different coefficients for Hq and 
H„i terms. One can, actually have a minimum gap that is larger than Eq. 6 by 
choosing a{s) = b{s) = \J s{s + 1) in the definition of the driving Hamiltonian 
Hd. In that case, the running time approaches 1 for the large N limit. 
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4 Conclusion 



In summary, we have presented the resuhs of an analysis of a generalized 
adiabatic quantum search algorithm. Our analysis shows that the optimal 
speed associated with Grover's search algorithm can be improved by a careful 
choice of the time-dependant Hamiltonian. This analysis does not deal with the 
complexity of implementing the driving Hamiltonian. Hamiltonian of adiabatic 
algorithms includes searched state \m) terms. The construction of these terms 
are thought to be analogous to having access to an oracle in original Grover 
algorithm. It seems that the time-complexity of searching in present approach 
is shifted from searching process to construction of the Hamiltonian. 
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